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Abstract 

The Misner initial value solution for two momentarily stationary black holes 
has been the focus of much numerical study. We report here analytic re- 
sults for an astrophysically similar initial solution, that of Brill and Lindquist 
(BL). Results are given from perturbation theory for initially close holes and 
are compared with available numerical results. A comparison is made of the 
radiation generated from the BL and the Misner initial values, and the phys- 
ical meaning is discussed. 



I. INTRODUCTION 

Black hole collisions are presently of great recent interest as one of the "grand challenges" 
in high performance computing The results of those studies, in turn, can be important 
to the understanding of detectable sources of gravitational waves [If. 

To the present date the only case that has been extensively studied is the head-on 
collision, from rest, of two holes starting with the initial value solution given by Misner 0. 
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The spacetime growing out of those initial data has been computed by the techniques of 
numerical relativity and has been studied by analytic means f^]- 

The initial value solution of Brill and Lindquist P) (hereafter "BL"), like the Misner 
solution, represents two initially stationary nonspinning holes. Neither solution contains 
any initial radiation of short wavelength compared to the characteristic size of the throats. 
Outside the horizon the two initial value solutions can be thought of as differing in the initial 
distortion of each throat caused by the presence of the other throat. There is in fact no 
solution of the initial value equations of general relativity that is uniquely singled out as 
representing two initially stationary holes. The Misner solution and BL solution are special 
only in their mathematical convenience, and in the topological properties of the geometry 
of the initial surface extended inside the throats. Specifically, the Misner solution may be 
thought of as having a two-sheeted topology. The two throats representing the two black 
holes connect an upper "physical" sheet to a single lower sheet isometric to the upper one. 
In contrast, in the three-sheeted BL solution, each of the throats connects from the upper 
sheet to a separate lower sheet. The isometry between the two sheets in the Misner solution 
results from an infinite series of image terms in the solution to the hamiltonian constraint. 
It is reasonable to expect that these terms might lead to additional gravitational radiation, 
not present in the BL solution. Other physical consequences of the image terms have been 
studied in Ref. 0. 

Here we extend the analytic study of collisions of holes to the case of BL initial data. 
There are two main justifications for doing this. The first is that analytic answers are a 
useful aid to development of the codes used in numerical relativity. The values reported 
here for radiated energy can be tested against numerical codes for evolution of axisymmetric 
initial data. For initially close black holes, it will be interesting to see whether those codes 
agree with the analytic answers as well as they do in the case of Misner initial data. 

The second reason for some interest in the evolution of BL data is the general question of 
the relationship of initial data to the generation of gravitational radiation. In astrophysically 
realistic problems the initial data will necessarily come from some approximation scheme, 
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such as post-Newtonian solutions. Such an approach is justified if the gravitational wave 
signal generated depends only on certain general features of the initial data and is insensitive 
to many details (e.g., topology). The comparison of the evolution of BL and Misner data 
gives us a simple model for studying this question, and an interestingly simple (though 
limited) answer. 

In the next section we describe the application of close-limit perturbation theory to the 
evolution of BL initial data. In Sec. Ill results are given for the radiation predicted by 
perturbation theory. These results are compared with available fully numerical results for 
the BL case, and are compared with analogous results previously reported for collision from 
Misner initial data. 



II. CLOSE-LIMIT PERTURBATION THEORY FOR BL INITIAL DATA 



Like the Misner solution, the BL geometry is conformally flat and takes the form ds^ — 
ds\^, where ds\^ is the line element for flat 3-dimensional space, and where $ satisfles the 

Laplace equation in the flat space. In terms of spherical coordinates i?, 6, (f), for ds^^, the 

Misner or BL metrics can be written. 



ds^ = ^\R,e;i^o) (dR^ + R^ 



de' + sin^ e 



2 n J ±2 



(1) 



For the BL geometry, the form of $, aside from a factor of 2, corresponds to the potential 
of Newtonian theory, with points of mass m at positions z — -^zq on the z axis: 



^BL = 1 + ^ 



m 



+ 



m 



(2) 



2 \ ^i?2 sin^ e + {Rcose- zof R^ sin^ 6 + {RcosO + zof 
For R> zq the square roots can be expanded in a power series in zo/R and the BL 3-geometry 



written as 



dsl^ = 



^-"^=0,2,...^-"^ 



(dR^^R^ [d^2^sin2^(i02]) , 
where the Pa are the Legendre polynomials, and where M = 2m. 



(3) 
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We next make a transformation of the radial coordinate i? to a new coordinate r, as 
if we were transforming, in the Schwarzschild spacetime, from isotropic coordinates to 
Schwarzschild coordinates: 



R = (v^+ Vr - 2M)V4 



(4) 



It is convenient now to rewrite the hne element for the 3-geometry as 
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M/{2R) 
1 + M/{2R) ^^^^ \MJ \R 



dr 



1 - 2M/r 



de"^ + sin^ e d(t? 

(5) 



where the meaning of R is given by 

The geometry in (^) reduces to the Schwarzschild geometry if the summation in the 
leading factor on the right is ignored. That summation, then, contains the information 
about the deviations from sphericity and is the starting point for close-limit nonspherical 
perturbation calculations [1^. In particular, the parameter e = ZQ/M can be considered an 
expansion parameter for perturbation theory. If, for each multipole index £, we keep only 
the leading order in e, the approximation to the BL initial geometry takes the form 



'BL 



dsv- ^ 



1 + 



2M/R 



1 + M/{2R) „J71 \MJ \RJ 



=2,4,... 



COS I 



dr"^ 



1 - 2M/r 



d9^ + sm^ 9 i 



(6) 



In principle, for each multipole index i, one can read off the metric perturbations (which 
are purely even parity) from (^, can construct Moncrief's |TT[] gauge invariant perturbation 
wave function tppert, and can evolve that wave function with the Zerilli equation ||T2||. In 
practice, this need not be explicitly carried out. There is a striking similarity between the 
expressions in (H)-® and the equivalent expressions for the Misner geometry p [§. The 
single difference is the coefficients in the series appearing in For the Misner initial 

geometry the coefficients are k^^Hq). The dimensionless quantity /iq parametrizes the initial 
separation of the throats, and the k's are functions given in Ref . Q . The single change 
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{z^/MY 4/€,(/xo) (7) 

converts to their equivalent form for the Misner case. This means, for a given £, 

that ippcvt for the BL case has precisely the same form as for the Misner case; the outgoing 
gravitational waves, according to perturbation theory, are identical in shape. They differ 
only in a multiplicative factor. Since power carried by outgoing waves is proportional to the 
square of ^pert the results for BL infall, for each £, can be found by multiplying the Misner 
results by [{z^/ / . We note in passing that the "forced linearization" procedure 
discussed in Ref. is, of course, also applicable to the BL data. 

This Misner-BL equivalence applies for any separation of the holes. For large separa- 
tions of the throats it is not surprising that the gravitational waves generated by BL and 
by Misner initial data should be similar. For small initial separations, however, there is a 
significant difference between the 3-geometries of Misner and BL, and it does seem strange 
that the gravitational waveforms should be identical. Furthermore, it is for close initial sep- 
aration that perturbation theory is most applicable, so the prediction of identical linearized 
waveforms is also a prediction about the actual waveforms. How can such different initial 
conditions give rise to identical outgoing waveforms? 

It is important to realize that the linearized outgoing waves are identical in form for each 
but the ratio of multipole contributions differs for BL and Misner. In Fig. 1 this difference 
in multipoles is shown quantitatively. For a given value of /io in the Misner geometry, an 
equivalent configuration for the BL geometry is defined by setting the quadrupole amplitudes 
of ippert equal, i.e., by setting {zq/M) = 2\Jk2{Jm))- The ratios of the BL amplitude to the 
Misner amplitude are then computed for £ = 4 and £ = 6. (These amplitude ratios are 
in fact simply 4K2(/io)^/^4(/^o) and 16K2(/io)^//t4(Ato).) At large separation the amplitude 
ratios approach unity; this shows that in the limit of large separation the external fields 
become identical in the two initial geometries. For small separations, however, the BL 
solution has a relatively smaller contribution due to higher multipole moments; its geometry 
is more quadrupole dominated. Though this is a relatively important difference in the initial 
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geometry near the throats, it is of httle importance for the gravitational radiation. Even for 
the Misner initial conditions, the radiation is heavily quadrupole dominated. It is possible 
that the lesson of this example has a broader generality: the outgoing radiation can be 
insensitive to many details of the initial data and even for strong field sources a knowledge 
of the quadrupole moment may be all that is needed. 

It is worth asking whether there is any deep physical meaning in the fact that the only 
difference between the BL and Misner linear perturbations is the ratio of the multipole 
amplitudes. This follows from the fact that for a conformally flat 3-metric, with the form 
(0), the factor $ satisfies the flat space Laplacian. If the solution is axisymmetric and 
asymptotically fiat it must be of the form ^{ae,/ R^^^)Pi{co^9)\ solutions can differ only 
in the values of the constants «£. So the striking similarity of the BL and the Misner 
perturbations is a direct result of the choice of the conformally fiat form (|l]). This choice is 
dictated by convenience, and need not be made in principle. For more general momentarily 
stationary initial geometries the linearized waveforms for each multipole will have different 
appearance. For example, one could generate valid initial data representing a Schwarzschild 
spacetime with a nonconformally-fiat perturbation by choosing an arbitrary (small) metric 
perturbation and solving the linearized hamiltonian constraint for the conformally-fiat part 
of the perturbation. The gauge-invariant function would then be computed from the full 
perturbation. 

III. RADIATION ENERGY: BL VS. MISNER 

The first, and most difficult, step in comparing radiation from the two initial value sets 
is to decide on the basis for comparison: How does one compare a BL problem with a 
particular value of Zq/M with a Misner problem of a particular /ig? At large separations 
it is not difficult; one can compare BL and Misner configurations in which the masses and 
separation of the holes are identical. For small separations, however, the separation of the 
holes is somewhat ambiguous. To deal with small, as well as large, separations we choose a 
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reasonably natural and convenient specific measure of the separation L: the proper distance 
along the symmetry axis, between the outermost disjoint marginally outer-trapped surfaces 
around each throat. (For zq/M less than about 0.4, a single apparent horizon encompasses 
both holes). The locations of the marginally outer-trapped surfaces was found using a 
standard shooting technique applicable to axisymmetric spatial slices 1113. We characterize 



both BL and Misner configurations with L/M, where M is the mass of the spacetime. It 
is, of course, interesting not only to compare the linearized predictions for BL against those 
for Misner, but also to compare both against the results of numerical solutions of the fully 
nonlinear field equations. For the Misner initial geometry the numerical results are known 
from the work reported in Ref. [HI. For BL initial conditions two data points are available: 



cases c2 and c4 of from Ref. [|14|. These numerically generated spacetimes have euclidean 
spatial topology, with initial data consisting of spherical (in the conformal space) collisionless 
matter configurations. When the initial configurations are sufficiently compact, the matter 
is all inside disjoint apparent horizons and the external 3-geometry is identical to the BL 
data. 

For clarity, the results are presented in three separate figures. Figure 2 shows the com- 
parison of perturbation results and numerical results for the Misner case. The perturbation 
energies {E/M ^ 0.025lK2(/io)) are those of Ref. f^, except that energy has been plot- 
ted as a function of L/M, rather than of fiQ. The numerical data are those of Refs. 0] 
and [^. Figure 3 shows the analogous results for the BL case, for which the energy is 
E/M ^ 0.025l[{zo/My/Af. The two "numerical" data points here are those of Ref. [|14|. 



Figure 2 shows that for the Misner case, linearized predictions begin to diverge from 
the fully numerical results at around L/M = 4. It is fortunate that the numerical results 
available for the BL case are for L/M in the range 3-4. From Fig. 3 we can infer that for 
L/M less than around 3, the agreement between linearized and numerical results is very 
good for BL collisions, and for L/M above 4 there is significant disagreement. In this sense 
there is little difference between Misner and BL cases. Figure 4 shows the perturbation 
theory comparison of Misner and BL cases. This figure shows that there is little difference 



between the predicted radiation when L/M is greater than around 2. It is, therefore, not 
surprising that the agreement between numerical and perturbation results, which breaks 
down well above L/M = 2, does not distinguish between BL and Misner collisions. It is 
also not surprising that in BL collisions, as in Misner collisions |0, the radiation is always 
quadrupole dominated. (The large values of hexdecapole energy in Fig. 3 occur only at 
separations large enough that linearized theory wildly overestimates radiation.) 

The results in Fig. 4 would seem to suggest that, for black holes initially close, BL 
initial conditions lead to less radiation than Misner black holes, as expected by the presence 
of image terms in the Misner solution. An alternative interpretation is that, for equal 
radiation, the initial separation of the apparent horizons is greater in the BL case than in 
the Misner case. Since equal radiation implies equal quadrupole moments, this means that 
the different multipole structure of the BL and Misner geometries makes the proper distance 
between apparent horizons larger in the BL case when quadrupole moments are equal. In 
this sense then. Fig. 4 is more of a depiction of proper distances than of radiation. 

This motivates asking whether there is a way of comparing BL and Misner scenarios that 
is better, or at least different, from using L/M. Another physically meaningful measure of 
how close the initial throats are is the gravitational binding energy. The gravitational 
binding energy is the difference between the ADM energy of an initial data set representing 
two black holes at finite separation and the energy of an initial data set with the holes 
infinitely separated (the sum of the bare masses of the holes). For BL data this is given by 



^ = --. (8) 

M 8zo ^ ^ 



For Misner data one has ||T5| 



Eb ^ En=iin - l)cschra/io 

M ~ E^=iCschn/io ^ ^ 

Radiated energy is plotted against binding energy in Fig. 5, but the results give a picture 
very much like that of Fig. 4. In particular, for small initial separations (tightly bound initial 
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configurations) there is less energy radiated from a BL collision than from a Misner collision. 
For large initial separations (small binding energies) the difference in radiated energy is 
small for configurations with the same binding energy. The BL and Misner cases become 
significantly different (say by a factor of 2) for binding energy (binding energy/M ^ —1.5) 
that corresponds roughly to the point {L/M 1.3) at which the BL and Misner energies 
separate in Fig. 4. 
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FIGURES 

FIG. 1. Ratio of amplitudes of ^pert for BL and Misner geometries. For equal amplitudes of 

£ = 2, amplitude ratios are shown for i = 4 and i = 6. 

FIG. 2. Gravitational radiation energy emitted during the head-on collision of two black holes 
starting from the Misner initial conditions. Results are shown for close-limit perturbation theory 
(continuous curve) and for numerical relativity (isolated points). 

FIG. 3. Gravitational radiation emitted during the head-on collision of two black holes starting 
from BL initial conditions. Results are shown for close-limit perturbation theory (continuous 
curve), and two values are shown from numerical relativity. 

FIG. 4. Comparison of perturbation theory predictions for radiated energy from Misner and 
from BL initial conditions. Results are given for both i = 2 and 4 multipoles. 

FIG. 5. Gravitational radiation emitted from Misner and BL initial conditions are plotted as 
a function of the binding energy of the initial configuration divided by the mass of the spacetime. 
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